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BONDED SOLIDS AT LOW TEMPERATURES 

Paul Whitney Sparks 
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Iowa State University of Science and Technology 

The linear thermal expansions of Si, Ge, GaAs, GaSb, InAs, and InSb 

were measured between 2°K and 30°K using a linear variable transformer 

technique . Sensitivities in length changes of 10 cm long samples of close 
o 

to 0.02A were observed. In general, the data show a region of positive 
thermal expansion at very low temperatures and this is followed by a region 
of relatively large negative thermal expansion at higher temperatures. 

From previous work it is known that at still higher temperatures the thermal 
expansion again becomes positive. The present data on Ge and Si agree 
satisfactorily with previous results. 

On a reduced plot, the features of the curves for the various solids, 
with the exception of Si, show a systematic dependence on 9 q . The meaning 
of these correlations and their lack of application to Si is not obvious 
from a theoretical standpoint. The values of 7 , the low temperature limit 
of the Gruneisen constant, determined for Ge and Si (0.66 and 0.44 respec- 
tively) do not agree with those calculated from the pressure dependence of 
the elastic constants (0.50 and 0.25). This is well outside expected errors. 
For T/9 q < 0.02, the thermal expansion coefficients of all the samples 

USAEC Report IS-T-114. This work was performed under contract 
W-7405-eng-82 with the Atomic Energy Commission. 
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vary as T . This is surprising since the specific heat does not show 

3 

T temperature dependence for Ge or Si except at much lower temperatures 
(T/0 < Q .01). As a result, y shows an initial decrease as the tempera 

ture increases. 
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INTRODUCTION 



A solid consists of a system of particles in stable equilibrium which 
vibrate constantly with small amplitudes. The free energy, 

F - U-TS, (1) 

of such a system is given by 

F - -kT InZ, (2) 

where Z is the partition function. In such a system the lattice potential 
energy may be written as 

3N 

0 - 1/2 S A q.q. + 0 q (V) , (3) 



where the q's are vector coordinates and the A.^'s are constants. Hence, 
the Hamiltonian may be written as 



3N ? 3N 

H - l/2m 2 Pj + 1/2 2 A q.q +0f(V) , (k) 

>1 J i,j J J 



where the P.'s are the momentum vectors and 0 is the lattice cohesive 
J o 

energy. This has a quantum mechanical solution with eigenvalues of the 
form 



e. * (n + 1/2) hv. + e Q , (5) 

where n is any given integer, v. is the frequency of the i th mode, and £ q 
is related to 0 q (V). Therefore, the partition function. 



Z - E exp (-e./kT) 

I 1 






( 6 ) 



( 



c 



2 






may be wr i tten as 
3N 

Z * S exp (-€ /kT) exp (-l/2hv./kT) E exp (-nhv./kT) 
i-1 ° ' n-0 ' 



J - 1 
=■ it exp (-E /kT) [1-exp (-hv./kT)] , (7) 

i-1 ° ' 



where E is a function of volume and contains the zero point and cohesive 
o 

lattice energies. Hence, 



3N 

F - E (V) + kT E ln[ 1 -exp (-hv./kT)] . (8) 



It is seen that this expression is of the form 

F ~ E o (V) + F*(V,T) , (9) 

where 

F*(V,T) =- kT E. f. (v./T) . (10) 

The quas i -harmon i c approximation which involves the assumption that the 
Vj's are a function of volume will be used here. From the definition of 
the free energy given above (equation 1) it is seen that 

F* * U*-TS , (11) 



and it is easily shown that 



U* 



, c>F*/T v 
Wt j \j 



K S 



df. 

d(v./T) V i * 



02 ) 



Similarly, the pressure, P, is given by 



. ( i£) 

WT 

SE -v r t 

- ( — -) - (^0 
V 'T v SV J 7 



03 ) 



A thermal pressure, P*, can be defined as 



p* 



-(— ) 
k av ' t 

df . 

-kT E (5^ ) 

I 



-k S 



df. dv. 
1 



t d(v./T) dV 



(14) 



Equations (12) and (14) can be used to show that 

df. 



w 

U. - k 



i ' N dTvviT v i J 



( 15 ) 



and 



,, d f . dv , 

a I I 

P i “ “ k d(v./T) d\T 



( 16 ) 



Hence, 



* 



it 

U. dv. 
1 1 



i v . dV 



U. d In v. 

1 1 



V d In V 



U" 

V 1 Y, 



( 17 ) 



where 



d 1 n v . 



The total pressure on the solid is given by 
P - P o (V) + P* (V,T) 

- P 0 (V) +S P* (V,T) . 

i 



09 ) 



The volume thermal expansion coefficient, P, and the isothermal bulk 
modulus, B.J., are related by the identity 

^ b t Wv * 

However, from equations (17) and (19) 



(££) * t— : ) 

wv 'st } \i 



y. 9U. 

v — L ( L\ 

. V 'dT ; V 

i 



y. C . 

' I VI 



( 21 ) 



t h 

where C . is the specific heat of the i mode. Therefore, 



P - 2 



Y-C . 

‘ I VI 



, V 



( 22 ) 



or 




( 23 ) 



where the Gruneisen parameter, y, is defined as 



f 



c 



» 
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7 



E C .y- 
v i ' i 



E C . 



v i 



(24) 



Equation ( 23 ) is relevant only to the lattice contribution to the 
thermal expansion coefficient; however, if there are other contributions, 
F^, to the free energy which are additive, one may write 

F * 2 F* (V,T) + 2 F (V) . (25) 

k k 



From equation (13) it is seen that 



(— ) 

^arv 



(S^L.) 
V avar v,t 



( 26 ) 



and from equation ( 20 ) 



P B 



T 



'3V9r V,T 



(27) 



Therefore, individual contributions to the thermal expansion coefficient 
may be defined as 




S ( -) 

“ \^1/^ T' 



T k 



5Vd rv,T 



(28) 



where B^. is the total bulk modulus at the volume and temperature, V and T. 

In the present experiment it is believed that contributions other than 
those due to the thermal energy of the lattice are of little importance and 
the most emphasis will be placed on equation ( 23 ). 

The Debye theory gives perhaps the most successful first order method 
for treating the thermal properties of solids. In this theory one considers 
a set of harmonic oscillators as in equation ( 8 ) and then makes an 



( 



( 






approximation to sum over the 3N modes. The approximation is that there is 
no dispersion between the wave vectors and frequencies in the solid., i.e. 
the velocity of sound is constant in the crystal. There is also a maximum 
frequency, v^, such that no frequencies higher than this value occur in the 
solid. The Debye theory then predicts that the specific heat of a solid 
will be given by 



hv 



T 3 e ° /T - 

c = 9R (5 f f — 

V \ H ( e x -l) 2 



4 , 

x dx 



(29) 



m 



where 0^ - — “ and is known as the Debye temperature. This approximation 
becomes very good at very low temperatures (T<0 q/ 5O) since in this range 
there is no dispersion in a real solid. In this case the Debye theory 
predicts that 



C v = 234R (T/0 O ) 3 . (30) 

A fit between experimental data and theory is usually made by considering 
the Debye temperature to be a function of temperature and in many cases 
specific heat data are presented as a plot of 0^ versus temperature. The 
limiting low temperature value of 0^ is usually designated as 0 Q . 

From the definition of y (equation 24) it is seen that at high tempera- 
tures (T»0 ^)y becomes constant since all modes have the classical value 
of the specific heat, C . = k, and 

VI 



Y = 



Y; 

? 3N 



(31) 



In the low temperature limit (T < 0 Q / 50 )y again becomes constant since only 
continuum lattice vibrational states are excited and there is no change in 
the relative weighting factors of the various modes. It is seen from 
equation ( 30 ) that the C^'s are proportional to r where the Y^'s are 
summed over the transverse and longitudinal acoustic modes so 



X) 



*> 



■ y = 



\ T ' 

E A. T 3 
k 



£v\ 

£ A, 



= const. , 



( 32 ) 



where the various A^ 1 s and s are constants. This low temperature limit 
may be expressed as 



= - 



d ln0 
d 1 nV 



( 33 ) 



where 9 q is the value of the Debye temperature at T = 0 K and may be related 
to the pressure dependence of the elastic constants. There is no reason for 
these two limits of y to be equal so the experimental y * s * n general a func- 
tion of temperature. For a typical solid such as NaCl the values of the high 
and low temperature limits of y are 1.51 and 0.93 respectively (l). 

At low temperatures , (9^/10), V , B^., and y are fairly constant so from 
equations ( 23 ) and ( 30 ) it is seen that 



(3 = const. T" 



(34) 



The thermal expansion coefficient of NaCI, for instance, is of the order of 
-10 3 

10 T . The present experiment, however, involves the thermal expansion 
coefficients of some diamond structure materials, for which the values of (3 
will be almost two orders of magnitude smaller than this. This arises 
because of smaller y 1 s and relatively larger 9‘s of these materials. It is 
therefore necessary to detect a relative change in the size of a solid of 
the order of a few parts in IQ 11 . 

1 SL 

Experimentally, the linear thermal expansion coefficient, o: = — (-^y) , 

o P 

is measured and for a cubic material Q! is related to 3 by 



a = P/3 



( 35 ) 



Experimental methods based on x-ray diffraction or optical interferometry 
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techniques do not have the sensitivity necessary to measure expansion 

coefficients as small as 10 ^/°K. Three very different types of apparatus 

have been developed which do have this sensitivity. First, a double grid 

arrangement has been used with one grid moving with respect to the other 

(2). The relative transmission of light through this combination is 

measured using a differential photocell, and results are quoted to a 
o 

fraction of an A. A modification of the above principle which uses a 

tilting mirror to shift the optical image of one grid with respect to the 

other has been tried at room tempera tures. The detection of sample move- 
_ qo 

ment of 10 A has been claimed ( 3 ), but this method has not been adapted 
to low temperature measurements. 

A second method involves plating the top surface of a sample and using 
this surface as one plate of a variable capacitor in a bridge circuit (4). 
This method is one of the most sensitive now in use for measuring low 
temperature thermal expansion, and changes in length of a few tenths 
angstrom can be detected. The major limitation of this method involves 
the stability of the capacitor plates and the need for small plate separa- 
tion. 

The method used in the present work consists of using a linear 
variable transformer to measure the expansion of the sample (5). As the 
sample changes length the motion is transferred to the secondary coil of a 
linear variable transformer and is measured by the use of a mutual induc- 
tance bridge. Changes in length can be detected with a precision of 
o o 

0.02A and a mean scatter of jt0.02A is maintained where necessary over a 
range of one angstrom. 

Prior to 1958 negative thermal expansions were observed in such solids 



( 



( 
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> 



as silver iodide (6), silicon (7), vitreous silica (8) , zinc blende (9), and 
a-uranium (10). This behavior in these materials is significant in that 
the negative thermal expansion is a result of their crystal structure and 
is strictly a lattice phenomenon. This is unlike magnetic materials such 
as invar where an additional magnetic term is added to the free energy to 
produce a net negative expansion. In 1958 Gibbons ( 1 1 ) measured the thermal 
expansion of Ge, Si, and InSb and found negative expansions for all of these 
materials at low temperatures. Following a suggestion due to Blackman (12) 
Gibbons concluded that negative y. 1 s for these materials could result 
since all of these materials have fourfold coordination (tetrahedral ly 
bonded), loosely packed structures, and are strongly covalent. 

These measurements were extended by Novikova over the next few years 

in an extensive investigation of the thermal expansion of diamond structure 

materials from room temperatures to as low as 20°K. Si ( 1 3) , Ge (14), and 

diamond (15) were measured first and the results confirmed those of 

Gibbons. The i soelectron i c series, Ol-Sn, InSb, and CdTe ( 1 6) ; and GaAs, 

ZnSe (17) and Ge (l4) were measured and it was found that the temperature 

at which o: becomes negative increases as the value of a . becomes smaller 

min 

when passing from a group IV material to a group II -VI material. This led 
the author to conclude that "... this result corresponds to the increase of 
the strength of the ionic components of the bonds." (17). In 19&3 
Novikova measured the thermal expansions of AlSb, GaSb, ZnTe, and HgTe (18) 
with the same general results. 

Carr, McCammon, and White (19) performed some very fine measurements 
on Ge and Si at lower temperatures (6°K to 100°K) with 0.5 angstrom sensi- 
tivity. This work demonstrated the presence of a region of positive 



( 
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expansion in these two materials from 0°K up to T/0 ^ 0.02$. The exi s- 

o 

tence of this low temperature positive region had been suggested by 

Daniels (20) using equation (32) and taking a weighted average of the 

pressure derivatives of the elastic constants based upon the acoustic work 

of McSkimmin (21). The thermal expansion data were not sufficiently 

precise, however, to give the details of the positive region although the 

negative expansion region was investigated in some detail. 

The most recent thermal expansion experimental work completed on 

diamond structure materials was again done by Novikova in a study of the 

effect of adding impurities to Ge (22) and Si (23). It was found that 

1 7 3 

electronic impurities with a concentration greater than about 10 /cm 
would increase by a few percent the absolute magnitude of the thermal 
expansion coefficient and the transition temperature where o: became nega- 
tive. The same behavior was shown to be given by a dislocation density of 
7 3 

the order of 10 /cm . This further confirms some work done on active 
impurities in Ge by Keyes (24). 

The theory of negative thermal expansions in general and that of 
diamond structures specifically has been studied for several years and 
although considerable insight has been achieved over this period there is 
unfortunately no complete theory to describe the details of this pheno- 
menon. Blackman (12) demonstrated that a negative thermal expansion 
coefficient was feasible in both the sodium chloride and zinc blende 
structures. From this work the phenomenon of negative thermal expansion 
passed from being an abnormal behavior of materials to possibly the normal 
behavior for materials with a loosely packed structure. 

As pointed out by Daniels (25) the existence of a region of negative 



< 



< 



thermal expansion is not the only anomaly of these materials. At high 
temperatures the value of OC i s almost an order of magnitude lower than 
expected of most materials. Also the specific heats of these materials 
exhibit large minima in the value of 9^ at about T * 0^/5. 

A more sophisticated model must be used to fit the existing data and 
Bienenstock (26) used a nearest-nei ghbor plus a dipole-dipole interaction 
model to make a one parameter fit to the experimental y's of Ge. A 
successful fit was achieved., but it is difficult to assign any physical 
significance to the one parameter. 

The purpose of this experiment is to investigate the low temperature 
thermal expansions of some diamond structure solids primarily in the low 
temperature region where their expansions are positive. It was hoped to 
determine if these materials behaved at very low temperatures as Debye 
solids (CC * AT ) and if any correlations and explanations can be given 
for their behavior in the temperature region where their expansions become 



anomalous 
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DISCUSSION OF APPARATUS 

The apparatus in general is based on that described by Carr (5) with 
several major modifications which will be discussed below. Figure (l) is 
a schematic representation of the primary components of the equipment and 
their relationship to each other. This equipment may be divided into four 
basic systems: sample suspension, sample heating, sample expansion measure- 

ment, and thermometry. These systems must respectively support and isolate 
the sample from its ambient surroundings, vary the temperature of the 
entire sample uniformly above the ambient temperature, measure the change 
in length of the sample, and measure the temperature of the sample. 

Polished sapphire disks are glued to both ends of the sample and to 
the base of the sample chamber. The disks at the bottom of the sample and 
at the base of the sample chamber are indented and a small sapphire ball is 
placed between them to provide a gimbal bearing. This bearing provides for 
the transverse support of the bottom of the sample and also corrects for 
small mistakes in alignment. Transverse support for the top of the sample 
is accomplished by attaching three nylon threads about one-fourth of the 
way down from the top of the sample, passing them over small pulleys, and 
attaching them to small springs. 

At the top of the sample a fourth sapphire disk is glued to a copper 
strip which leads to the bath. When the sample is mounted this disk is in 
contact with the disk at the top of the sample and a quartz rod is placed 
between the copper strip, i.e. the top of the sample, and the secondary coil 
of the linear variable transformer expansion sensing device. Since there is 
a vacuum in the sample chamber, the polished surfaces of the disks provide 
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Figure 1. General electrical and mechanical 
schematic drawing of the thermal 
expansion measuring apparatus 
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for a high thermal impedance and at the same time maintain good mechanical 
contact throughout the expansion system. 

Sample heating is accomplished by passing a small current through a 
manganin wire heater which is wound on and glued to the sample. This heat 
will raise the temperature of the sample uniformly (if the thermal impedance 
of the sapphire disk junctions is much greater than that of the sample) 
until an equilibrium is reached where the heat delivered into the sample is 
equal to the heat loss through the ends of the sample into the bath. 

A linear differential transformer is used to measure the sample 
expansion. The secondary coil is in mechanical contact with the sample by 
means of a quartz rod. The primary and secondary coils are wound in such 
a way that when an alternating current is passed through the primary coil 
a voltage is induced in the secondary coil which for small displacements 
from the center position is proportional to the change in position. This 
induced voltage is sensed as a change in the coupling between the primary 
and secondary coils and is measured using a mutual inductance bridge. 

The temperature of the sample is measured with an Au-Fe versus Cu 
differential thermocouple (27). One junction of the thermocouple is 
attached to the sample and the other junction is in the helium bath. 

Since the bath temperature can be determined from vapor pressure readings 
(it is held at a constant temperature below the helium X point with a 
large Cartesian manostat) the absolute temperature of the sample can be 
found directly from the thermocouple voltage,, which is measured by a 
potent iometric method. A detailed discussion of the various equipment 
components will be given in the following sections. 



( 



( 










f) 



16 



Sample Chamber and Sensing Coil 

Once the approximate size of the sample has been decided upon, four 
basic problems must be solved: a method for transmitting the motion of 

the sample to the secondary coil, a method for supporting the secondary 
coil, a method for supporting and insulating the sample from bath temper- 
ature., and a method for introducing the electrical leads into the sample 
chamber. 

In this experiment the secondary coil is placed inside of the sample 
chamber and a direct mechanical linkage is provided between the sample and 
the secondary coil (see Figure 2). This method eliminates the need for a 
diaphragm as used by Carr (5). The use of a diaphragm allows the secon- 
dary coil to be placed directly in the helium bath so that it is maintained 
at a constant temperature. However, if a sample has a large thermal 
expansion, it will tend to pull away from the diaphragm upon cooling, 
while a sample with a small thermal expansion will tend to rupture the 
diaphragm upon cooling. The present method offers greater flexibility, at 
the expense of requiring rigorous thermal anchoring of the secondary. 

The secondary coil is supported from a phenolic shield inside the 
sample chamber by three nylon threads glued to- the top and bottom of the 
secondary coil and to the shield. The secondary coil is wound without a 
coil form and a quartz tube with one end sealed and indented is fixed 
axially within the coil. The quartz tube is glued with G.E. 7031 adhesive 
to the midpoint of the secondary coil since this minimizes the change in 
the balance point of the mutual inductance bridge upon cooling from room 
temperature to liquid helium temperatures. This also minimizes any effect 
due to extraneous heating. A 2mm quartz rod is glued to the sample thermal 
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Figure 2, Cutaway drawing of the sample chamber 
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ground (copper strip) which is in direct mechanical contact with the sample 
The other end of the rod is placed in the indentation in the quartz tube so 
that the secondary coil is connected to the sample. 

This type of construction eliminates the use of a diaphragm and the 
associated loading problems, since the secondary coil is inside of the 
sample chamber. It is necessary, however, to insure that the temperature 
of the secondary coil stays extremely constant and close to the bath 
temperature, since heating of the secondary coil can affect the balance 
point through thermal expansion or susceptibility changes of the secondary 
coil. Since the secondary coil is in a vacuum in the sample chamber a 
mechanical heat conductance path to the bath is necessary to remove heat 
generated in the secondary coil by eddy currents. This is accomplished by 
placing three 0.002 in. copper strips longitudinally between the windings 
of the secondary coil, soldering them to another copper strip that leads 
down the phenolic shield in the sample chamber, and attaching this strip 
to a 0.06 in. copper rod which leads to the bath. To prevent heat from 
the sample from affecting the secondary coil, the sapphire disk at the top 
of the sample and the base of the quartz rod are attached to a sample 

thermal ground which consists of a 0.003 in. copper strip soldered to 

another copper rod which leads directly to the bath. In addition, a 0.01 
in. copper wire is glued to the quartz spacer and to one of the copper 
rods. To check these precautions Fe-Au versus Cu differential thermo- 
couples were temporarily installed on the secondary coil and on the sample 
thermal ground referenced to the base of the sample chamber. The tests 

were conducted with a primary current of 100 ma and a sample heater input 

of 20 milliwatts (sufficient to raise the sample temperature to about 30°K) 
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The temperature of the secondary coil rose 0.6°i< due to eddy currents from 

the primary field and was unaffected by the sample heat. The temperature 

of the sample thermal ground rose 2°K due to the maximum sample heating 

power and was unaffected by the primary current. Calculations show that 

these values of heating of the secondary coil and sample thermal ground 

should produce effects that are two orders of magnitude less than the 
o 

least count (*0.01 A) of the experiment. 

As briefly described above the sample is supported and insulated 
through the use of sapphire disks and nylon threads. An important require- 
ment in this experiment is that the thermal conduction of the sample be 
adequate to prevent thermal gradients from developing across the sample. 

If the thermal impedance provided by the sapphire disks is high, then less 
power can be used to elevate the temperature of the sample to the desired 
level and, therefore, there will be a smaller temperature gradient across 
the sample. In this experiment this condition is checked by comparing 
data taken with the complete heater in use with data taken with only one 
half of the heater in use (a center tap on the heater is provided for this 
purpose). if there are important thermal gradients in the sample these 
two sets of data will not agree. 

A polished sapphire disk 1mm thick and 1/4 in. in diameter is glued 
to each end of the sample with G.E. 7031 adhesive. Another set of disks 
is glued to the base of the sample chamber and to the bottom of the sample 
thermal ground. If two hard polished surfaces are in contact in a vacuum 
a large thermal impedance is developed across the junction (28). There- 
fore these disks provide for a good thermal break and at the same time 
maintain mechanical continuity from the base of the sample chamber, through 



( 



( 



21 



> 



the sample, along the quartz spacer, and to the secondary coil. The disks 
at the bottom of the sample and at the base of the sample chamber are 
indented and a 1mm sapphire ball is placed between them forming a gimbal 
bearing. This bearing provides transverse support for the base of the 
sample and also adjusts for small misalignments. 

Three size A nylon threads are tied to the sample about one-quarter 
sample length from the top of the sample to provide transverse support for 
the top of the sample. These threads lead over small pulleys and are tied 
to springs constructed of spring bronze. The springs are attached to 
traveling screws so that after the threads are tied small adjustments may 
be made in the alignment of the sample. The pulleys and the traveling 
screws are glued and soldered respectively to a copper shield which 
surrounds the sample. This shield, which is made from 99.999% pure copper, 
is used both for sample support and for mounting the secondary phenolic 
shiel d. 

With the incorporation of a thermocouple as the primary thermometer 
in the system, it became necessary to use a continuous-lead vacuum-tight 
connector to insure that thermoelectric voltages are not developed where 
the electrical leads leave the sample chamber. Hollow pin connectors, 
which were tried previously, were found to develop leaks a majority of the 
times that they were cooled to liquid nitrogen temperatures. 

A connector was developed for this purpose (see Figure 2) by fitting 
a 1/8 in. stainless steel tube into a hole drilled in a 1/2 in. section of 
3/8 in. diameter nylon rod. These two pieces are glued together with 
Armstrong A-4 epoxy adhesive after which a 1/32 in. hole is drilled in the 
nylon through which the wires are passed and glued. It is very important 
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that all surfaces including all of the wires are thoroughly cleaned and 
completely wetted with the glue before the joint is made. After the glue 
joints have been cured the wires are passed through a copper tube into the 
sample chamber and the stainless steel tube is soldered to a copper 
adapter with Wood's metal. 

The remainder of the sample chamber consists of the glass envelope 
(1.75 in. O.D.) which has a 1.5 in. copper to glass feathered seal on one 
end and a piece of optical glass blown into the other end. This envelope 
is Wood's metalled to the base and the sample chamber is now vacuum-tight. 
The primary coil is placed over the outside of the glass envelope., is held 
transversely with two phosphor bronze springs., and is held longitudinally 
with a phenolic “cap". A screw which is in this cap pushes on the flat 
top of the glass envelope., and it is used to make gross adjustments (to 
about +0. 1mm) in the balance point. 

Since the sample lengths varied from 7 cm to 13«5 cm it was necessary 
to build two sample chambers. The first sample chamber was used for 
sample lengths of the order of 10 cm (although one sample was only 8.2 cm 
long) and was designed as described above and as shown in Figure 2 . For 
the longer and shorter samples a second sample chamber was constructed 
along the same design with the exception that the sample shield was 
lengthened to 5"3/8 in. and small hooks were placed inside the shield 1-3/^ 
in. from the bottom. For the longer samples the loading procedure was the 
same as described above. However., for the shorter samples the nylon threads 
were passed over the pulleys and down through the hooks on the inside of 
the sample shield and then connected to the sample. This allowed the sample 
to be supported at a point 3"l/2 cm from the bottom of the sample chamber. 
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A slit was cut in the sample shield near the top of the sample for the 
sample thermal ground. Once the short sample was mounted a long quartz 
spacer was installed and the remainder of the assembly was carried out as 
described above. 



Dewar System 

For an experiment of the present sensitivity dewar design becomes 
quite important for the elimination of magnetic., eddy current, and vibra- 
tional effects. The introduction of a glass dewar system eliminated the 
magnetic and eddy current problems and a simple innertube suspension system 
has eliminated a greater part of the vibrations (see Figure 3). 

The present glass dewar system was designed with a large bulge to 

contain sufficient liquid helium for a ten or twelve hour run. The long 

tail of the helium dewar was designed so that very long samples (up to 

20 cm) could be used if desired to increase sensitivity. In this case 

there still would be sufficient liquid helium above the sensor coil for a 

five or six hour run. It was also decided to seal the nitrogen dewar so 

that the liquid nitrogen could be pumped down to its freezing point. The 

bubbles from liquid nitrogen at one atmosphere cause vibrations of about 
o 

5A. If the nitrogen is pumped down to its freezing point and sealed-off, 
no boiling will occur until its vapor pressure has again reached one 
atmosphere. It generally takes five hours after pump-down is completed to 
boil again. 

The dewar system and liquid coolants combined weigh about 200 pounds. 
This system is suspended rigidly from a 1/2 in. thick aluminum plate which 
is supported on two small innertubes which lie on the support frame. These 
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Figure 3 
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Schematic drawing of glass dewars and 
dewar system support 
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innertubes act as very effective dampers for frequencies of the order of ID 
to 20 c/s which form the major vibrations external to this experiment. The 
pressure in the innertubes is of the order of five psig and may be varied 
to adjust the damping constant of the system. 'All pumping and sensing 
lines have at least a small section of rubber hose in them to maintain the 
flexibility of the system. 



Heater Supply 

The temperature of the sample is raised by passing a small current 
through a 40 Q heater wire (number 44 manganin) which is wound around the 
sample and is attached with G.E. 7031 adhesive. The power used is of the 
order of 20 milliwatts at 30°K. A center tap and each end of the heater 
wire are connected to a low current power supply in such a manner that 
current can be passed through either the whole heater or only through one- 
half. This is done so that thermal gradients across the sample can be 
detected by comparing "end only" measurements with "total" measurements. 

The current supply is a transistorized current source which will 
provide direct currents from 0 to 400 ma in four ranges to a low resistance 
load. The lower ranges are used in taking data at liquid helium tempera- 
tures while the upper ranges are used for systems checks at room and liquid 
nitrogen temperatures. A change in line voltage of 20% or a doubling of 
the load will change the current by on 0.1%. This stability is necessary 
since small fluctuations in the heater current will produce significant 
fluctuations in temperature. 
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Mutual Inductance Circuit 

The variable transformer consists of a fixed primary coil, a movable 
secondary coil, and an additional pickup coil wound around the primary 
turns (see Figure 4) so that a modified mutual inductance bridge can be 
used. The primary coil is wound in three sections where the ends (1/4 of 
the coil length each) are wound in one direction and the center half is 
wound in the other. This provides a field configuration as shown in 
Figure 5. The important features of this configuration are that the field 
approaches zero rapidly at the ends of the coil thereby reducing eddy 
currents in the sample and sample chamber which may produce heating or 
erroneous signals in the mutual inductance bridge. The middle and end 
“humps" are flat to about 0 . 1 % over about 1mm. The pickup coil is wound 
directly over the primary coil in the same manner to provide maximum 
1 inkage. 

The secondary coil is made in two halves wound in opposite directions. 
The length of the secondary coil is such that when the coil is centered 
the ends of the secondary coil lie on the end “humps" of the primary field. 
Therefore, when the secondary is situated midway with respect to the 
primary field there is no net flux linkage. As the secondary coil is moved 
with respect to the primary coil a small distance a voltage is induced in 
the secondary coil whose phase differs by 1 80° depending upon the direction 
of movement. The magnitude of this voltage is proportional to the coil 
off-set. Since the sample is attached mechanically to the secondary coil, 
the voltage produced will be a direct measurement of the sample expansion. 

The axial magnetic field at the center of the primary coil divided by 
the primary current can be expressed by 
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Linear variable transformer and winding 
data 
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Figure 5. Primary field configuration with magnitude 
corresponding to a primary current of 
100 ma 
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where ji. o is the permeability of free space., is the number of turns per 
unit length of the prima ry coil, and y is a geometrical factor (about 0.7 
for the present coil). B q for the present coil is about 300 gauss/ampere. 
The mutual inductance between the secondary and primary coils is 



A n 

M = (2z'B - 2z'B.) , 

l v o r J 



( 37 ) 



where A is the mean area of the secondary coil, n is the number of turns 
s s 

per unit length of the secondary coil, is the magnetic field at the end 
of the secondary coil (i.e. at the “hump"), and z‘ is a small displacement 
from the null position. However, B q -Bj can be expressed as FB q , where r 
is another geometrical factor, so 



M(z') = (2/1 )(A n B Fz 1 ) = 2A n n ji, 2 



s s o 



s s p r o 



( 33 ) 
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since Yf = 1 . This coupling is about 8 x 10 henries/A for the present 
transformer. The induced voltage in the secondary coil then is 



E(z') = -M = -2coA n n ll ! z 1 
dt s s p*<D p 



( 33 ) 



where a) is the angular frequency of the primary current. Thus the sensi- 
tivity is given by 



— = 2a>A n n jj. I 

d z 1 s s p*o p 



(40) 



For the present experiment where U)/2it is 2$0 c/s, I is 100 ma, and A g , n g 

and n are given in Figure 4, 
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It was apparent in the earlier work (5) that the stability and range 

of the various "standard" mutual inductances of the conventional Hartshorn 

bridge caused serious limitations to both the sensitivity and the precision 

of the measurements. For these reasons a ratio transformer was incorporated 

into a modified mutual inductance bridge (see Figure 6). In this way the 

mutual inductance between the primary and secondary coils is cancelled by 

the fraction of the mutual inductance between the primary and pickup coils 

which is determined by the ratio transformer setting. The pickup coil and 

the ratio transformer therefore replace the variable inductance of the 

1 0 

Hartshorn bridge and allow a relative range of 10 which is accurate to 

0.001% of the full scale reading, which does not drift with temperature, 

and which can not be changed by vibration or jolting. 

For ease of calibration, it is convenient to have an effective bridge 

length of more than one millimeter. At the same time the least count of 

o g 

the experiment is 0.01A so a range of at least 10 is necessary. It was 

1 0 

decided to design for a range of 10 with a full scale value of one centi- 
meter to insure that at least the one millimeter range was achieved. The 
ratio transformer used in the present experiment is a modification of two 
Gertsch model 101 1R seven decade ratio transformers connected to provide a 
range of 10 (see Figure 6). The inputs of these two transformers from 
the pickup coil are connected in parallel. The high side of one of the 
transformers is the high side of the ratio transformer complex. The output 
of the other transformer is reduced by 1000:1 transformer handwound around 
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Figure 6. Mutual inductance bridge circuit 
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a Supermalloy toroid core. The high side of the output of the toroid 

transformer is the low side of the ratio transformer complex. In this way 

7 -3 

the second ratio transformer has a range or 10 which starts at 10 of the 
first ratio transformer, so a total range of 10^ is achieved. Since the 
mutual inductance between the primary and the pickup coils should equal the 
mutual inductance between the primary and the secondary coils when z 1 is 
lcm, the number of turns on the pickup coil is 



N 



2A n z* 
s s 



pu ~ A 



(42) 



pu 



or about 1000 turns. 

The alternating voltage source is a variable frequency power oscillator, 
Optimation model AC-15. The primary loop is series-tuned so that a lower 
output voltage can be used from the oscillator for a given primary current. 
The gain and the frequency of the oscillator are stable to within 0.01% 
over a twenty-four hour period. This is necessary since a fluctuation in 
either of these settings will affect the balance point significantly due 
to capacitive pickup between the various coils and to other non- i dea 1 ized 
elements in the circuit. Two dropping resistors (0.1 to 10 ohm) are 
included in the circuit to cancel resistive components. The amount of 
voltage across these resistors is regulated by a 1KH rheostat which controls 
the current in the resistive loop. 

The null voltage in the secondary loop is detected by a narrow-band, 
low noise, lock-in amplifier-detector (29). The amplifier has a low noise 
input (+ 0.05 vrms with a one second time constant) and a low microphonic 
response. The least count signal to the amplifier is increased by a factor 
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of 500 through the use of an S.I.E. model 0155^9 input transformer. The 
null is observed on a d.c. ammeter output from the lock-in detector. An 
oscilloscope is used to monitor the phase of the signal. 

In working with a mutual inductance system where changes of a few 
parts in lo' ' of the flux linkage are important, many problems are 
encountered which would not be significant in a system of lower sensitivity. 
In this experiment, considerable time was spent in isolating and elimina- 
ting second and third order effects. It is worthwhile to discuss some of 
these problems and their solutions even through in many cases it is not 
known which of several possible causes led to spurious effects. These 
problems will be discussed in two categories: magnetic and eddy current 

effects and grounding effects. 

In an ideal bridge, a change in the primary current should not cause 
a change in the balance. However, if a magnetic material is present in the 
vicinity of the sensor coil this is not the case since the permeability 
ji = ji.(H) . For example, in the first stages of this experiment a 1% change 
in primary current caused a shift in the balance point of 10 ^ on the 
ratio transformer (full scale being unity). Subsequently it was found that 
ni chrome wire had mistakenly been used instead of a non-magnetic wire to 
suspend the sample chamber from the top of the dewar and that these wires 
ran along the outside of the sensor coil. When these wires were replaced 
with linen thread the current effect was reduced by a factor of 500. This 
current effect was reduced even further by replacing the Kovar seal on the 
glass envelope by a copper- to-gl ass "housekeeper" seal. Since the oscil- 
lator now in use has an amplitude stability of 0.01% the current effect 
(which appeared as drifts or noise) is no longer a factor in the 
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measurements. 

A more subtle problem was encountered when the required sensitivity 
and noise level were obtained, A 10 cm long sapphire sample was installed 
and an experimental run was attempted to compare the results with data 
taken two years ago on an interim apparatus. The results obtained indicated 
that there was a large anomalous "expansion" as shown in Figure 7° The two 
most significant points of interest in these results are the regions of 
apparent negative expansion at about 1 1 » 5°K. and 2S°K. The lower region was 
found to be primary field dependent., but not a function of bath temperature. 
The upper region, on the other hand, was found to be independent of primary 
field, a function of bath temperature, and larger than the lower effect by 
a factor of 3- 

Au-Fe thermocouples had been installed to monitor the temperatures of 
the secondary coil and the sample thermal ground. Part of the anomalous 
behavior may have been a result of the localized magnetic moments in the 
thermocouples so these were removed. Also, stainless steel screws which 
were used to adjust the tension on the nylon sample support threads were 
replaced with brass and the sample thermal ground was simplified. When 
these changes ware completed another run was made and the upper dip had 
disappeared, the lower dip was slightly reduced, and with the exception of 
the lower region the data from this run and that of the previous experiment 
were identical. 

No further adjustments were successful in eliminating the lower dip so 
a new sample holder was constructed. Another run was made with results 
which were within 1% of those obtained two years earlier. It is believed 
that this success was due to the use of Cd-3i non-superconduct i ng solder for 
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Figure 7. Anomalous results due to various magnetic 
a. id eddy current effects in the thermal 
expansion apparatus 
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all connections in the vicinity of the sensor coil and to the use of high 
purity copper (39.999%) in all places where copper was used in the sample 
chamber. This can be important as it has been found that small amounts of 
magnetic impurities in copper can result in a rapidly varying susceptibility 
as a function of temperature (30). It is not certain if these two possi- 
bilities were the cause of the lower dip., but they seem to be most probable,, 

In a mutual inductance circuit such as the one described here one must 
be very careful with the grounding system. In principle it should be 
possible to establish the system ground at any point in the circuit as 
long as the rest of the system is isolated properly. However, due to the 
juxtaposition of large grounded metal objects such as the dewar frame and 
the electronic rack and their capaci ti ve coup 1 ing with the elements of the 
circuit, it was found that the system ground must be placed at the low side 
of the detector input. If this was not done large spurious signals 
resulted in the circuit. 

It is also necessary to eliminate all high empedance paths between 
any high potential element and ground or another element. It was found, 
for example, that the five-way connectors on the inputs and outputs of the 
ratio transformers had a capacitance of 300ji.f between them. These were 
removed and replaced by BNC connectors. All lead wires were thoroughly 
shielded. It also is necessary to eliminate even minor ground loops. 

These wiring faults would cause anything from complete signal saturation to 
changing the ratios on the ratio ttansformers. 
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Thermometry 

This experiment was designed to operate in the temperature range 
between 2°K and 30°K» In order to obtain an accuracy comparable to that of 
the calibration of the mutual inductance bridge (0.2%) it is necessary to 
have an absolute temperature determination of +0.01°K„ In this experiment 
this criterion is met by the use of an Au-Fe (0.07 atomic percent) versus 
Cu differential thermocouple referred to the bath temperature ( 27 ). 

This thermocouple has a sensitivity of about l4y,v/°l< throughout the operating 
temperature range so voltages of up to 450p.v must be measured to +0. lp.v. 

These voltages are on the lower ranges of commercial potentiometers, and 
because of the small absolute voltages involved., it was decided to construct 
a direct current nulling circuit (see Figure 8) , 

One junction of the thermocouple is glued to the sample and the other 
junction is in the helium II bath. The two copper leads of the thermo- 
couple lead to room temperature and the thermocouple voltages would normally 
be detected between them. However, the nulling circuit is constructed by 
placing a small resistor, r^ (0.6443 + 0 . 0003 ^), in series with one of the 
two copper leads that are then connected to a Kiethley model IpO microvolt- 
meter. An independent pair of leads is used to pass a small current 
through this resistor to produce sufficient voltage to n bucl<-out n the 
thermoelectric voltage. The null is detected on the microvoltmeter. A 

standard resistance, R r - 1000Q, is placed in series with the current 

s 

source and the voltage drop across this resistor,, (as measured by a 
Keithley model 662 differentia] voltmeter),, is proportional to the thermo- 
couple voltage,, V^,, and is given by 
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Figure 3 0 Thermocouple measuring circuit 



THERMOCOUPLE CIRCUIT 








The ratio r d /R s is determined at 1.9°K by replacing the null detector 
with a microvolt potentiometer, and then comparing the zero AT voltage 
drops across the thermocouple leads (or r^) and for a fixed current 

through r^ and R s » 

It is seen that will be a maximum of about 0„7 volts (when the 
sample temperature is 30 °K ) and must be determined to + 155 p.v for +0.01°K. 
Since all other connections in the circuit with the exception of the sample 
junction and the connection at the null indicator are in liquid helium II 
only small (less than 1.0 pLv) constant thermoelectric voltages are developed. 
The thermoelectric voltage at the connection at the null indicator is kept 
small since this connection is thermally anchored to a copper plate and 
encased in styrofoam. The total stray voltages indicated at the microvolt- 
meter are of the order of lp.v., and do not vary more than +0. ]p„v over ten 
hours. These voltages are measured when the system is in thermal equili- 
brium with the bath (no sample heat) and used to correct the differential 
thermocouple voltages. 

The current source must be stable to 0.01% over short periods of time 
(during a measurement). Since the maximum current is only 0.7 ma, a bank 
of nine mercury cells connected in series is quite successful in providing 
the voltage for a standard d.c. current supply. 
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CALIBRATION 



Since the total linear bridge length was designed to be well in 
excess of one millimeter and calculations indicated that the primary field 
should be linear to better than 0.5% over this range a direct calibration 
was used. This method involves using a traveling microscope slide to move 
the secondary coil with respect to the primary coil (see Figure 9). This 
instrument can be read with care to +1 micron. The sensing coils are placed 
in liquid nitrogen to take into consideration the thermal contraction of 
the coils between 300 and 77°Kj and to increase the sensitivity of the 
bridge by reducing the resistance of the coils. One calibration run was 
performed at liquid helium temperatures as a check and the same results were 
obta i ned. 

For the “top half" of the coil, the calibration of the sensing coils is 
AL = 1 . 130(+0.00l) An cm where n < 1 is the ratio transformer reading, i.e. 
1.130 cm for the full bridge range. An = 1 (see Figure 10). The “top half" 
of the coil is defined such that the center of the secondary coil is above 
the center of the primary coil. The calibration of the “Dottom half" is a 
few percent less than the “top half", but since it is not as linear (about 
0.4%) the “bottom half" is not normally used to take data. 

It should be emphasized that the use of this calibration depends upon 
the fact that the variable transformer remains linear (or becomes more 
linear) as the range of motion decreases. Therefore, the calibration 
obtained over a larger range can be trusted when making small measurements. 
This can be checked somewhat by running a sample at different bridge settings 
and checking the results. This was done and no discrepancies were found. 
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Figure 9. The calibration apparatus 
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EXPERIMENTAL PROCEDURE 

The various components of the system have been described above. The 
i nter-rel at ionsh i p of these can be appreciated better through a description 
of the procedure involved in the performance of an experiment. Three steps 
are involved; sample preparation and loading, systems check and helium 
transfer, and the method of taking data. 

The samples used in this experiment vary in length from 7 cm to 13»5 
cm and are about 1/4 in. to 3/8 in. in diameter. As described above the 
10 cm samples are used in sample chamber 1 while the shorter and longer 
samples are used in sample chamber II. In preparation, the sample ends are 
polished normal to the sample axis and the sapphire disks are glued to the 
ends with G.E. 7031 adhesive. The heater wire is wound astatically on the 
sample and glued, and then the thermocouple is attached in the same way 
making sure that good thermal contact is maintained along two or three 
inches of the sample. All electrical connections are made, the sample is 
mounted on the sapphire ball, the quartz spacer rod is installed, and the 
sample chamber is assembled. 

After the primary coil has been installed and the pumping line attached, 
the sample chamber is pumped and flushed with nitrogen gas several times 
and then leak-detected. A rough balance is set with the primary coil and 
the whole apparatus then is placed in the dewar where it rests on the 
bottom. This is done so the measurements will not be affected by small 
changes in the dewar length. 

Once the apparatus has come to equilibrium at 77°K and before liquid 
helium is transferred the heater, mutual inductance, and thermometry 
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circuits are checked out once again* It usually is not necessary to 
rebalance the primary since there is sufficient range on the rough balance 
at room temperature to compensate for any differential expansion (which 
could be~0«>3mm). Once these prechecks are completed the helium dewar, 
vacuum chamber,, and the sample chamber are sealed-off. Liquid helium then 
is transferred and after another set of systems checks are made, the bath 
is pumped down to 1*89°K and regulated by a cartesian manostat. After 
about one hour the whole system has settled down and data can be taken* 

In this experiment there are two distinct types of data taken; at the 
lowest temperatures a continuous record of AL versus temperature was taken 
every few tenths of a degree. At higher temperatures five points are 
taken over a small temperature range (about 0*3°K) and fitted with a straight 
line to obtain the thermal expansion coefficient, a, directly, where 



a = — (vy) (see Figure 11)* The CL points were taken every one or two 
degrees and AL, a length change, was taken between these points to give a 
continuous AL versus T plot if needed* This could be obtained by integra- 
ting the CL vs. T data also* Although these two types of data give different 
quantities the method of taking them is quite similar* 

First a base point is established by having the system come to equili- 
brium with a given heater current (at the lowest temperatures this would be 
= 0) * At this point the thermocouple voltage and the mutual inductance 
are recorded. A small current in excess of the base point current is then 
put into the heater, the temperature rises, and another point is recorded. 

It was found to be convenient to use the Kiethley 150 voltmeter readings 
for the small thermocouple difference readings as data are taken between 
base points, so as the heater current is increased from a base point the 
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thermocouple voltage difference is read on the Kiethley up to 10 + 0 o 1 
(about 0o7°K)o The linearity of the meter readings of the microvoltmeter 
used in this experiment was checked and found to be 1%. When the thermo- 
couple reaches 10 *j,v a new base point is established and the process 
continued,, Normally a reading is taken every 2 jjlv so that about seven 
points per degree are taken* To get the total length change the base point 
lengths then are added. In this experiment this method was continued until 
the thermal expansion of the sample became negative and then CC points were 
taken* The method for obtaining CC was the same as at lower temperature 
except that five points were taken above the base point at thermocouple 
voltage intervals of 1 p,v* 

In taking all of these points at least two sets of data were taken at 

every base point (up to twelve sets at the lowest temperatures) and the 

o 

data averaged. This allowed for a mean scatter of + 0*02A where needed 

o 

even though the least count of the apparatus was only + 0*05A* Also some 
data were taken with only one half of the heater in use to check the temp- 
erature distribution across the sample* Points also were taken both up and 
down in temperature to insure that there was no hysteresis in the expansion* 
In all cases the base point was returned to repeatedly insuring that there 
were no drifts or jumps in the balance point while data were being taken* 
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RESULTS 



In this experiment thermal expansion measurements were performed on 
S i , Ge, GaAs, GaSb, InAs, and InSb. These materials were chosen since they 
form a group centered around Ge on the periodic chart and it was hoped that 
correlations would be found which would assist in understanding the behavior 
of these tetrahedral ly bonded solids. The sample properties and sources are 



given in 


1 Table 1. The 


original 


20.5 cm long Ge 


sample was cut 


into two 


pieces., 


one approximately twice 


the length of the other, so that data from 


samples 


of d i f ferent 1 


engths could be used to check the apparatus for syste- 


matic errors. The samples were 


prepared and data were taken as described 


previously. The actua 


1 experimental data are p; 


-esented in the 


append I x* 


Table 1. 


Sample properties and 


sources 














Net Carrier* 






4°K Length 


D i ameter 


Res istivi ty 


Concentration 




Sample 


_ . (cm) 


(cm) 


(Q~cm(300 o K) ) 


(atoms/cm^) 


Source 


Si 


13«65 


~0.9 


150 


-,o 15 


Texas Inst. 


Ge 


13.46 


-1.2 


48 


<io 17 


Mono S i 1 i con 




6.99 


~1.2 


48 


<10 17 


Mono S i 1 i con 


G aAs 


8.26 


~0 0 6 


0.002 


-2 x 10 18 


Texas Inst. 


GaS b 


10.06 


'* s -0 0 6 


0.68 


~io 17 
-1 /* 


Texas Inst. 


1 nAs 


10.01 


'■"'Oo 6 


0.012 


-2 x 10 ° 


Texas Inst. 


InSb 


10.12 


-1.0 


0.024 


<1.6 x 10 16 


Ohio Semi- 












conductors 



'The various uncompensated impurity concentrations were determined by 
comparing the observed resistivities with those of doped samples. Data for 
Ge., Si, GaSb, and InSb were obtained from Hannay (3l)j while the information 
for GaAs and InAs were obtained from unpublished literature supplied by 
Texas Instruments, Inc. In addition, the reflectivities of Ge and Si were 
kindly measured by Professor D. W. Lynch, of this laboratory (private 
communication 1966). These results indicated only that impurity carrier con 
centrations were less than 5 x 10^7. Helicons were detected in a piece of 
the InSb sample by J. L. Stanford, also of this laboratory (private communi- 
cation 1966) indicating the stated purity. 
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As was expected from earlier work, as the temperature was increased 
each sample indicated an initial region of relatively small positive thermal 
expansion which was followed by a region of relatively large negative expan- 
sion (see Figures 12 to 17). In about the first two-thirds of these posi- 
tive expansion regions the length changes of the samples are proportional 
to T . This would be expected from a simple application of the Debye 
theory. At higher temperatures the region of negative thermal expansion 
continued to the highest temperatures reached in this experiment ( 30 °K - 
35°K) , although a minimum thermal expansion coefficient, a . , was observed 
in all cases except for Si. 

The mean scatter in the AL's in the low temperature positive thermal 

o 

expansion region was as low as *0.02A for InAs and InSb and for all samples 

4 

the T coefficients were determined with an accuracy of from 1.5% for Si to 
6% for InSb (see Figures 12 to 17). The smoothed values of a and Yy which 
are given in Table 2 are believed to be accurate to about 2% up to 22°K. 
Above this temperature the Au-Fe thermocouple calibration was obtained by 
extrapolating data taken at lower temperatures, so theAT's in this region 
could be in error by as much as ten percent. A proper calibration at these 
higher temperatures is planned. The results obtained for Ge and Si in this 
experiment compare favorably with those of Carr et al. (19) and a comparison 
of the data is shown in Figures 12 and 13 and in Table 2. It is seen that 
the data for Ge agree quite well except at higher temperatures where the 
thermometry in this experiment is in doubt. The data for Si do not agree 
as well, but the data of Carr et al. are in some doubt due to sample 
problems (19). It is difficult to compare the results obtained for GaAs, 
GaSb, and InSb with the data of Novikova (16, 17, 18), since that author's 
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Figure 12„ Linear thermal expansion of silicon 
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Figure 14. 



Linear thermal expansion of GaAs 



o< 



6j 



) 







AL, 



I0 3 & 



0 



4 



( 






64 



) 



) 













I 



l 



Figure 15* 



Linear thermal expansion of GaSb 
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Figure 16. Linear thermal expansion of I nAs 
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Figure 17. Linear thermal expansion of InSb 
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Table 2. Smoothed values of a and y 



S i • Ge 



r 

T 


a( 0 K _1 ) x 10 8 


y ! 

T T j 


CC( 0 K) _1 x 10 8 


■>v 


0 


0.0 


0.44 


0.0 


0.66 


2 


7.0 x lO -4 


0.44 


6.2 x 10" 3 


0.65 


4 


5.6 x 10~ 3 


0.44 


4.9 x 10" 2 


0.63 


6 


1.9 x 10 -2 


0.43 


1.7 x 10 -1 


0.60 


8 


4.5 x lO -2 


0.42 


0.40 ( 0.35) 


0.53 ( 0.45) 


10 


8.8 x 10 -2 


0.41 i 


0.50 ( 0.45) 


0.36 ( 0.25) 


12 


0.152 ( 0.15)“ 


0.38 ( 0.38) i 


0.45 ( 0.60) 


0.12 ( 0.15) 


14 


0.251 ( 0.12) 


0.34 ( 0.18) j 


0.21 ( 0.25) 


0.03 ( 0.035) 


16 


0.32 ( 0.10) 


0.26 ( 0.09) i 


-0.33 (-0.45) 


-0.03 (-0.035) 


18 


0. 18 (-0.08) 


0.05 (-0.04) 


-1.2 (-1.3) 


-0.06 (-0.065) 


20 


-0.18 (-0.45) 


-0.07 (-0.17) 

| 


-2.3 (-2.3) 


-0.03 (-0.080) 


22 


-0.70 (-0.85) 


-0.15 (-0.21) 


-3-7 (-3-6) 


-0.10 (-0.092) 


24 


-1.45 (-K45) 


-0.22 (-0.25) 

j 


-5.2 (-4.9) 


-0.11 (-0.10) 


26 


-2.6 (-2.6) 


-0.28 (-0.33) ! 


-6.8 (-5.6) 


-0.11 (-0.091) 


28 


-4.3 (-3-8) 


-0.33 (-0.36) , 


-8.2 (-6.1) 


-0.11 (-0.084) 


30 


-6.7 (-5.2) 


-0.37 (-0.38) 


-9.4 (-7.1) 


-0.11 (-0.084) 


i 


a = 0. 88 (+0.. I6)xl0 


" 12 T 3 ,T<15.9°K 


a « 7.80(+0.32 )x10~ 12 

j 


T 3 jT<9. 2°K 




G a As 




1 

| GaSb 




T 


|a(°K _ 1 ) x 10 8 
i 


Tt 


CO 

0 

X 

1 

0 

Is 


'Vt 


0 


0.0 


0.56 


0.0 


0.32 


2 


6.9 x 10“ 3 


0.56 


ro 

1 

0 

X 

-tf 

• 

CO 


0.31 


4 


5.5 x 10" 2 


0.54 


6.7 x 10" 2 


0.28 


6 


1.9 x 10 -1 


0.52 


2.26 x 10 -1 


0.24 


8 


4.4 x 10 _1 


0.43 


-0.18 


-0.10 


10 


0.45 


0.25 


-2.4 


-0.32 


12 


0. 10 


0.01 


-7. 1 


-0.43 


14 


-0.95 


-0.11 


-12.9 


-0.45 


16 


-2.7 

j 


-0.18 


-18.8 


-0.44 


18 


-5.1 


-0.22 


i- 


-0.43 



f 
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Table 2. Continued 





GaAs 






GaSb 






T 


a(°K _1 ) x lo 5 


^T 


aCV 1 ) 


CO 

2 

X 


^T 




20 


- 7° 7 


-0.24 


-23 




- 0.42 




22 


-10.5 


-0.26 


-34 




- 0.40 




24 


-14.0 


-0.27 


-38 




- O.38 




26 


-18.0 


-0.28 


-42 




- 0.35 




23 


-22.6 


-0.29 


-43 




- 0.32 




30 


-25.6 


-0.29 


-41 




- 0.23 




32 


-25.4 


-0.26 










34 


-21.8 


-0.21 












a = 8.40(+0. l6)xlo 


~ 12 T 3 ,T<8.0°K 


a = 10.7(±0. 1 6) xl 0 


" 12 T 3 ,T<6.6°K 






InAs 






1 nSb 






T 


oo 

0 

X 

1 

o 

a 


^T 


1 

O 

'a 


3 

x 10 


'Yt 




0 


0.0 


0o 33 


0.0 




0.19 




2 


1.1 x lo" 2 


0.31 


1.23 


-2 

x 10 


O 

O 

OO 




4 


0.0 


0.00 


0 

0 

0 

0 




0.00 




6 


-0.25 


-0. 17 


-1.00 




-0.46 




8 


-1.6 


-0.39 


-9.4 




-0.73 




10 


-5-3 


-0.43 


-32.1 




-0.86 




12 


-12.5 


-0.52 


-42.2 




-0.88 




14 


-22.7 


-0.55 


-63.7 




-0.83 




16 


-33-0 


-0.58 


-35.3 




-0.88 




18 


-45.2 


-0.60 


-105 




-0.83 




20 


-58.7 


-0.63 


-125 




-0.87 




22 


-68.5 


-0.65 


-144 




-0.55 




24 


-33.0 


-0.66 


-162 




-0.83 




26 


-97.5 


-0.68 


-175 




-0.81 




23 


-112 


-0.69 


-181 




-0.77 




30 


-127 


-0.69 


-182 




-0.72 




32 


-141 


-0.70 


-180 




-0. 66 




34 


-153 


-0.69 










36 


-151 


-0.68 












a = 13.6(±0.7)x10 


" 12 t 3 ,t<3.5°;< 


a - 16.4(±0.8)x10‘ 12 T 3 ,T<3.0°:< 






-.'The values in 


parenthesis are 


those of Carr et al. 
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limits of error are understandably large in the region of overlap (20°K - 
35°i<) and discrepancies as large as a factor of three occur. 



T able 3. Phy sical constan t s reduced to 0 D X 









Si 


Ge 


GaAs 


GaSb 


! nAs 


! nSb 


e 


°K 




645 a 


374 a 


34S b 


267° 


247 b 


206 b 


O 

c n 

C 12 

c 44 


io 11 


3 

dynes/cm 


16. 77 c 


1 3» IS 0 ' 


12.26 s 


8.97 f 


cn 

(VA 

O 

CO 


6 . 66 h 


io 11 


3 

dynes/cm 


6.50 c 


4.95 d 


5.71 e 


4. 12 f 


4.53 9 


3.35 h 


io 1 1 


3 

dynes/cm 


8. 04 c 


6.84 d 


6 . 00 e 


4. 51 f 


3.98 s 


3 . 1 4 h 


B I 


io 11 


3 

dynes/cm 


9.85° 


7.65 d 


7.89 s 


5.95 f 


6 . 02 s 


4.45 h 


a 


o 

A 


3 


5.44' 


5.67' 


5 . 66 s 


6. 1 ]■* 


6 . 08 s 


6 . 50 s 


P 


g/cm 


2. 33 1 


5.32 1 


5-31 1 


5 . 62 1 


5.87' 


5 . 78 ' 


V 


3 

cm 




12.06 


13.59 


13.55 


17.13 


16.74 


20.45 


m 

AC, ,/AP 
AC 12 /4P 
AC 1 ,/AP 


amu 




28.03 

4.23° 

4.20 c 

0.75 c 


72.59 
5 . 12 k 
4.45 k 
1 . 3 1 k 


72.32 


95.73 


94.87 


1 18.29 



a Flubacher et al. (32) 



Verma et al . (33) 
c McSkimmin et al. (34) 
d McSkimmin (35) 
e Garland et al . ( 38 ) 

f 

Einspruch et al. ( 37 ) 

S Ger 1 i ch (33) 
b Potter (39) 

Aigra i n (40) 

•* Giesecke et al. (41) 
k Kopperman et al. (42) 

The results of thermal expansion experiments are expressed best in 
terms of the Gruneisen constant, y, which is defined by equation 23, 



30® T V 




( 23 ) 
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ihe low temperature value 7 q can be obtained from the limiting region where 

-3 

both the thermal expansion coefficient and the specific heat vary as i . 

The specific heats were calculated using equation 30 and 0 q values obtained 
from the elastic constants of the various solids. The data necessary to 

calculate these 7 Q ‘5 are given in Table 3° In general., y is proportional 

3 3 3 3 

to 0:0 /T and the errors in both Q!/T and 0 are believed to be the same 
o o 

order of magnitude. 

In order to determine the temperature dependence of the Gruneisen 
constant (7^.) it is necessary to know the temperature dependence of the 
specific heats. Neither of the other parameters (V , B^) exhibit a signifi- 
cant temperature dependence. Unfortunately^ low temperature specific heat 

data exist only for Ge and Si (32) (where T/9 > 0.01) and for InSb (43) 

o 

(T/0 > 0.03). The specific heats of the other materials (44) have been 
o 

measured only down to 14°1< (T/9 q > 0*05). These data are not adequate for 
these 7^ calculations. An extrapolation of the existing specific heat data 
was performed using a reduced plot of the existing data of the form 9/3 

o 

versus T/9^ (see Figure 18). The extrapolation of the incomplete data was 

made by comparison with the reduced plots for Si and Ge using the values 

in Table 3- One striking feature (shown in Figure 18) is that at no 

3 

temperature does 0 appear to be constant so the true T region for the 

specific heat appears only well below T/9 = 0*01* It is possible that 

o 

errors in 9 as large as five percent may result from this procedure., but 
the variation in 0 is more than 30% so at least a realistic value of 9 (and 
therefore C^) may be obtained by this method* The actual values of were 



o Reduced 0 versus reduced T curves used for 
o 

the extrapolation of data to low temper- 
atures with 0 as given in Table 3® 
o 3 



Figure 18 
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taken from tabulated values of the Debye function. Smoothed values of y.. 
are given in Table 2 , while a reduced plot of the y' s (y/y Q versus T/8^) 
is given in Figure 19. 
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DISCUSSION OF RESULTS 

The results of this experiment confirm the existence of a low temper- 
ature region of positive thermal expansion for the materials measured. The 
general behavior for all these materials is similar to that shown for Ge in 
Figure 20. A relatively small positive thermal expansion at very low 
temperatures is followed by a region of relatively large negative thermal 
expansion and then by a region of positive thermal expansion for T/6 q > 0.2. 

The thermal expansion coefficient, Ct 3 has been shown to be proportional 
to T J at the lowest temperatures (i/Q o < 0*015) as would be expected from 

the Debye model* However ^ the heat capacities of neither Ge nor Si obey a 
3 - 

T law even down to i/G < CL 01 and this factor has been used in the extra- 

o 

polation on Figure 18, Hence., must increase constantly toward y q as the 

temperature decreases from T/G < 0.015. The deviation of heat capacity 

3 

from a T behavior is only a few percent for T/0 < 0*01., so the variation 

in for Si below 6°K is less than the experimental error (about 2%) * 
However., a definite decrease in can be observed in all cases before the 

3 

(X = OL r region has been passed (see Figure 19). The variation in shape 
o 

with S q for the various reduced y versus T curves can be seen in Figure 19. 

There appears to be a definite correlation with the obvious exception of Si. 

Sheard (45), Daniels (20), and Collins (46) all have predicted 
el 

*s for Ge and Si using the pressure dependence of the elastic constants 

cl cl 

for these materials. For Ge, y Q is 0.50 and for Si, y Q is 0.25, the 
calculation being based on the low temperature elastic constant vs. pressure 
data of McSkimmin (21) which is shown in Table 3« These calculations 
disagree with the values for Ge (y - 0.66) and Si (y Q = 0.44) obtained in 
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this experiment. At very low temperatures these materials should behave 
as continuum solids so the calculations should apply as well as for other 
solids (46). These discrepancies are well outside of the expected errors 
in both the high pressure elastic constant data and the low temperature 
thermal expansion data, so no reason can be given at present for these dis- 
crepancies. The data of Carr et al. (1.9) agree with the data of this 
experiment within the combined experimental errors (see Figures 12 and 13) 
making an error of 30 to 50 percent in the thermal expansion data unlikely. 

Another source of error in this experiment could involve an electronic 
or carrier contribution to the elastic constants as described by Keyes (24). 
Butj as can be seen in Table 1 the carrier concentrations of our samples 

] 7 9 

(with the exception of GaAs) are all 10 atoms/cnr or less and impurity 
effects should be negligible., 

In an attempt to determine the existence of correlations between the 
physical features of the thermal expansions of the various materials in 
this experiment, a tabulation of these features (as defined in Figure 21) 

is shown in Table 4, These features include T. and T ( , the two temperatures 

L n 

at which a = 0; T^, the temperature at which the expansion first deviates 

from a i ‘ behavior: and T . , the temperature at which the minimum thermal 

min ' 7 

expansion coefficient occured, Also are listed AL ^ , the height of the 

max 

low temperature positive thermal expansion normalized to a 10 cm sample 

lenqth; y . , the value of the minimum y„; y , the low temperature limit of 
3 'mi n ‘ > o 

y-, and Cl the coefficient of in the expression for CC ( =Cl T^) below T n . 
‘To o u 

It can be seen from Table 4 that with the exception of Si there 

appears to be a definite correlation between these physical features of the 

thermal expansions and the value of 0 . T . (which is roughly constant) is 

o min 



Figure 21« General behavior of the thermal expansion 
for a typical tetrahedra 1 1 y bonded solid 
with physical features denoted 
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Table 4 C Physical features of the thermal expansion s 







Si 


Ge 


GaAs 


GaSb 


1 nAs 


InSb 


T, 


°:< 


18.6 


14.7 


12.2 


7.7 


4.6 


3-9 


t h 


°!< 


120 a 


43 b 


55° 


52 d 




67® 


T o 


°K 


15.9 


9.2 


8.0 


6.6 


3-5 


3-0 


T . 
min 


°i< 

0 


75 f 


32 


31.2 


28 


35 


30 


AL 


A 


17.9 


32.5 


22.4 


6.65 


0.79 


0.60 


max .. 


a . 


10 V K 


-5 


-0.95 


-2.63 


-4.40 


-15.3 


-18.3 


min - 


\ii n 




-0.43 


-0. 1 1 


-0.29 


-0.45 


-0.73 


-0.89 


Y Q 




0.44 


0. 66 


0.56 


0.32 


0.33 


0.21 


a 

0 


T 3 x 10' 12 1°i< 


0.88 


7.80 


8.40 


10.7 


13.8 


16.4 



Nov i kova et al, ( 1 3) 



Novikova (14) 

C Novikova (17) 

^ Nov i kova et al. ( 1 8) 

6 Novikova (16) 

^ Carr et al, (19) 

an exception. Therefore., a plot of the Y 0 ' s and 5 was made as a func- 

tion of 9^ (see Figure 22). With the exception of Si both the Y Q , s and the 

y . l s appear to be a linear function of 0 „ If one considers these mater- 
r min rr o 

ials to be members of a “group 11 (excluding Si) then it may be possible to 
predict that any member of this “group 11 which has a value of 9 q below 
137 (+10) °K will not have a region of positive thermal expansion at very low 
temperatureso Also., any member which has a value of 9 above 405 (+10) °K will 
not have a region of negative thermal expansion. Hence., the value of 9 q for 
a member of this “grouo“ can be used to obtain both y and y . . Unfortuna- 
telyj at the present time there are no criteria for the conditions necessary 
to become a member of this “group 11 and the part which Si takes is unknown. 

It would be valuable to obtain the low temperature thermal expansions of such 
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materials as A1P(0 q -~ 5S0°i<) , A!As(0 ~ 420°K) , and CdTe(© o ~ 140°K) to 

investigate the “anomalous" behavior that Si takes in these correlations 

and to test the predictions made for a member of the "group". Thermal 

expansion measurements have been made on CdTe by Novikova (16), but these 

measurements do not extend to sufficiently low temperatures or are not 

sufficiently sensitive to test the suggestion that will be close to 

zero or even slightly negative in this case. It also would be valuable if 

more data existed for the pressure dependence of the elastic constants of 

other tetrahedral ly bonded solids. Perhaps the discrepancies involved 

between the y 1 s so obtained and those of thermal expansion measurements 
r o 

could be resolved. This type of measurement also might show the same type 
of correlation as has been indicated in this experiment. Also some informa- 
tion on the y' s may be gained from phonon-assisted tunneling under pressure 
as investigated by Payne (47) as this experiment measures y. at the band 
edges . 

The low temperature limit of the Gruneisen constant., y Q , is given by 
the volume dependence of 0 q (see equation 33). de Launay (48) has prepared 
tables for the evaluation of 0 q using the relationship 



e* - & <t > 3 (J r> 3/2 - V '<*•*) 

° 4«V k p 



(44) 



where f(s,t) is a relationship between the elastic constants, V is the 
molar volume, and p is the density. Using equation 33 j it follows that 

d 1 n0 

o 



d 1 nV 



89 



1 d lnC ^ _ 1 _ I 8 1 nf (s, 1) 

2 d InV 6 38 InV 



( 45 ) 



Daniels (20) computed for Ge and Si using this relationship and observed 
that the third terms were small and did not vary much between the two 
materials. Therefore, since the second term is constant the first term is 
of primary importance. This term may be expressed as 



d ln c ii4 i b t dC A4 



2 d In V 



2 C 44 8P 



( 48 ) 



Since the ratio for the materials measured in this experiment varies 

from 1.12 (for Ge) to 1.52 (for InAs), the large variation in y (a factor 
of three) must arise from the pressure dependence of C , ;l , . Unfortunately, 
it is difficult to develop an expression for dC^/8P from first principles 
and to relate this expression to 0 . At the present time there is no way 
to justify theoretically a linear relationship between the y's and 6 , or 
to show why Si differs from the other five materials. 
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CONCLUSIONS 

The results of this experiment have demonstrated the presence of a 
low temperature region of positive thermal expansion in the materials 
measured. It seems apparent that all tetrahedral ly bonded solids, except 
possibly those with a very low 9 q , should exhibit the same behavior. The 
experimental behavior of these solids is understood qualitatively in terms 
of the pressure dependence of the cu vs. k relationship for the lowest 
transverse acoustic mode. The pressure dependence of the elastic constants 
(for Ge and Si) have shown that 7 > 0 for small k while phonon-assisted 
tunneling experiments have shown that 7 < 0 at the band edge (47). Hence, 
the shape of this dispersion curve must change appreciably with pressure, 
while the other phonon branches show "normal" (7 > 0 ) behavior. 

The band structures of Ge and Si are fairly well understood in basic 

terms, but the pressure dependence of these bands involves too many 

unknowns and further theoretical work is necessary. Any theoretical basis 

for the relationships between 7 ‘s and 7 . 's and 9 is nonexistent at 

r o min o 

6 1 

the present time. Also, the discrepancies between the 7 0 's and 7 ^ ‘s for 
Ge and Si are not understood unless, perhaps, continuum theory is not 
applicable to the tetrahedral ly bonded structures even at low temperatures. 
The need for further work on the pressure dependences of the elastic 
constants of these materials is obvious. 

To study the thermal expansions of diamond structure materials in more 
detail would require more sensitivity. The present apparatus has been 
used close to its maximum sensitivity, and to increase the precision of the 
experiment it might be worthwhile to consider the use of superconducting 
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wire in the linear variable differential transformer. With the use of 
superconducting windings the primary current, and therefore the sensitivity 
could be increased by a factor of 100 as the present primary current is 
limited by the helium boil-off rate in the dewar. Also if the secondary 
coil were constructed with superconducting wire the matching into the 
detector-amp 1 i f i er would be almost unlimited,. Unfortunately, at the presen 
two problems exist. First the eddy-current heating produced by the primary 
current in the sample chamber could be uncontrollable. Second the use of 
low frequency alternating currents in superconducting wire is not fully 
understood, and hysteresis losses and locked-in flux are greater than can 



be tolerated. 
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APPENDIX 

Experimental Data 

These tables contain the data taken in this experiment. The column 

headings stand for the Au-Fe thermocouple reading (TC) , the temperature 

(T) obtained from our present ca 1 i bra t i on ; and the change in length of the 

sample (AL). These data were obtained by establishing first a base point 

which corresponds to a fixed value of the heater current ! u , and sample 

n 

temperature (initially this would be with 1 = 0) • At this point the 

n 

thermocouple, TC, and ratio transformer, n, readings would be recorded. 

The heater current then would be increased slightly (by switching to a 
different control on the heater current supply), the TC and An would be 
recorded, and then the sample would be returned to the base point. At low 
temperatures readings would be taken for AT's of approximately 0.1 or 0.2 
degrees. At intervals of one or two degrees a new base point would be 
established and the procedure repeated. While taking data between base 
points values of An would be taken with both increasing and decreasing 
temperature and system checks would be made as described previously. 
Therefore, all data points are taken as differences from the previous base 
point. All points were taken at least twice (at lower temperatures they 
were taken as many as twelve times) and these values were averaged and 
are shown in order of increasing temperature in these tables. At higher 
temperatures the same procedure was followed except that base points were 
taken every 2 degrees or so and the set of points were taken at 0.1 
degree intervals. 

The low temperature data are shown in these tables with the base points 
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underlined and the groups separated by a space. The base points of the 
high temperature data have a line both over and under the value and a 
space is also provided between groups. At lower temperatures theAL.'s 
between each point are summed giving a continuous AL as the temperature 
increases. At higher temperatures only the base AL.'s are summed since 
the AL.'s within a set are used only to determine a. In some instances 

i 

a base point at higher temperatures was not taken and in these cases the 
base AL is indicated by a dash. Since this was done many times for 1 nAs 

(the first sample measured) an additional set of data was taken, with 

o o 

only K89 K as a base point, up to 36 K. 
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